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An explicit solution is found for the most general independent correlation func- 
tions in lattice QCD2 with Wilson action. The large-A^ limit of these correlations 
O , may be used to reconstruct the eigenvalue distributions of Wilson loop operators 

for arbitrary loops. Properties of these spectral densities are discussed in the region 
/3 < /3c = i 
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^ . In view of the renewed interest in the master field approach to large-A^ matrix field 

theories [00, it may be convenient to improve our knowledge of those solvable models that 
^ \ will eventually become suitable laboratories for testing new methods and techniques. Lattice 
QCD2 with free boundary conditions is such a solvable system. In principle all correlation 
^ ■ functions can be computed by a proper gauge choice leading to a complete factorization of 
^ . the functional integral and by knowledge of the properties of invariant group integration. 
D ! The underlying single matrix problem has been solved a long time ago , and at first glance 
^ ' not much can be added to what is already known. 

However, as extensively discussed in Ref. [0], one of the problems that must be faced in 
^ ■ even simpler models (like a model of n independent Hermitian matrices) is that of finding 
a master field description of correlations involving different, and in general non-commuting, 
matrix field variables. The purpose of the present letter is that of establishing a benchmark 
for this problem, in the form of an explicit large- expression for such correlations at arbi- 
trary distance, and studying some properties of the eigenvalue distributions of the (master) 
Wilson loop operators associated with these correlations. 

Our starting point is the axial gauge formulation of lattice QCD2 with Wilson's action 
presented in Refs. [|,|],^. In this formulation the expectation value of a Wilson loop is 
reduced to that of an invariant product of field variables in a one-dimensional principal 
chiral model with free boundary conditions. Let's therefore focus on these models, defined 
by the lattice action 

S = 2f3Nj2^eTT[UiU}^^] (1) 

i 

where i are the lattice sites. By straightforward manipulations one may show that the most 
general nontrivial correlation one really needs to compute is 
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where / plays the role of the space distance and /c is a sort of "winding number" for the loop. 
In turn by the variable change Vi = Ui-iU} , the evaluation of Wi^k is reduced to computing 
the group integral appearing in the relationship 



JdV,...dVi j.TT{V,...Vi)'' exp(2iV/5EUReTrV,) 



[J dVexp{2NpReTTV)f 
This problem can be solved for arbitrary by a character expansion (see e.g. [^): 

exp[2iV/3ReTr\/] ex J] rf(,)%)(/?)xw(V^) , 
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where X{r) are the characters and d^r) the dimensions of the irreducible representations of 
U{N). Let's remind that the character coefficients 5(r)(/3) are explicitly known for all groups 
U{N). As a consequence of Eqs. we obtain 
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where we have introduced the coefficients Ck,(r) of the decomposition 

TrX'^ = J2Ck,(r)Xir)iX). 
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It is now an exercise in representation theory to show that, for all k < N, the only repre- 
sentations occurring in the r.h.s. of Eq. (M) are those in the form 



X(i,fc) = X(fc-i+l,p-l;0) ) 
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and one may show that 
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As a consequence we immediately obtain the final result holding for A; < A^ in all U{N) 
groups: 
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where we also know that 
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Further closed-form results may be obtained by restricting attention to large A^ and 
focusing on the strong coupling domain. When /5 < 1/2 and A^ is large B 
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Under the same assumptions therefore 
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and one can finally show that 
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We have explicitly checked that Eq. (|13D satisfies the Makeenko-Migdal equations for the 
large- limit of the corresponding Wilson loops in QCD2. 

Eq. ([T3|) is our main result, in that it is the quantity that a master field approach to 
QCD2 on the lattice should be able to reproduce for arbitrary / and k by substituting a 
single field configuration in the correlation functions. 

It is interesting to explore the properties of the Wilson loop operator at a distance I, 
and in particular its eigenvalue distribution, by constructing a generating function of the 
expectations Wi^k- 



Mx) = J2Wi,kx'' 
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From the Makeenko-Migdal equations it is possible to derive the Schwinger-Dyson equation 
satisfied by Ai{x) in the strong coupling regime (3 < 1/2 and in the large- iV limit 



{A{x)-l)Mxy-' = f3'x. 



(15) 



We checked explicitly the consistency of Eq. (|1^) with our solutions {^)- When / = 1 
one recognizes Gross- Witten's solution for the strong coupling phase of the single plaquette 
models. When / = 2 the solution 
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1 + AP^x + 1 
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is related to the generating function for the moments of the energy density 
^^^) = l + 2xfJ^ + 2 

k=l 



iV^^'i-/3x(K + KVi) 



^) = 1 + 2xp^ + 2Y,{(3x)^''W2,k = 2(3^x+ ^1+A(3^x^ . (17) 



Incidentally, Eq. (p!7| ) corrects a mistake in Ref. 0. 

More generally, for arbitrary /, we may reexpress our results in terms of the eigenvalue 
distribution pi{9) of the Wilson loop operator ni=i ^5 which is related to the generating 
function Ai[x) by: 
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and satisfing 



r depi{9) = 2Wi^, -1 = 1. (19) 

J —TT 

The study of the area dependence in the spectral density of continuum QCD2 was pi- 
onered a long time ago P,pD| and recently reconsidered by many authors. The lattice 
counterpart of this problem is the study of the dependence on I of the functions pi{6), which 
can be performed numerically, for arbitrary /, by summing up sufficiently many terms in the 
Fourier expansion of Eq. (pISl). In Figs. [| and || we plotted the spectral densities as functions 
of 6 for different values of / at fixed values of /3, and in particular when /3 = /3c. 

As one may notice, the property 

pK^) 1^ , (20) 

which can be independently derived from the confinement properties of the model |T^, is 
exhibited by our solutions. It is however quite interesting to notice that at /? = /3c and 9 = % 
the equation satisfied by the eigenvalue density is 

[l-2npi{n)][l + 2npi{n)f-' = 1. (21) 

Besides the strong coupling solution, which corresponds to the values plotted in Fig. |^, 
Eq. (0) admits a second solution pi^n) = (coinciding with the first one when / = 1,2). 
This phenomenon leaves open the possibility of a compactified solution in the weak coupling 
phase for arbitrary values of /. 
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FIGURES 

FIG. 1. Pi{0) at /3 = ^ and for various value of /. pi{—0) = pi{6) by symmetry. 
FIG. 2. Pi{0) at P = Pc = ^ and for various value of I. 
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